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Abstract 

We have shown by using the exact solutions for the two-electron system in a parabolic con- 
finement and a homogeneous magnetic field [l|, 0] that both exact densities (charge- and the 
paramagnetic current density) can be non-interacting V-representable (NIVR) only in a few 
special cases, or equivalently, that an exact Kohn-Sham (KS) system does not always exist. All 
those states at non-zero B can be NIVR, which are continuously connected to the singlet or 
triplet ground states at B = 0. In more detail, for singlets (total orbital angular momentum 
Ml is even) both densities can be NIVR if the vorticity 7(r) = V x ^j p (r)/n(r)^ of the exact 
solution vanishes. For Ml = this is trivially guaranteed because the paramagnetic current 
density vanishes. The vorticity based on the exact solutions for the higher \Ml\ does not vanish, 
in particular for small r. In the limit r — > this can even be shown analytically. For triplets (Ml 
is odd) and if we assume circular symmetry for the KS system (the same symmetry as the real 
system) then only the exact states with \Ml\ = 1 can be NIVR with KS states having angular 
momenta m\ = and 1 777-2! = 1- Without specification of the symmetry of the KS system the 
condition for NIVR is that the small- r-exponents of the KS states are and 1 . 

PACS numbers: 

31.15.E- Density-functional theory 

31.15.ee Hohcnbcrg-Kohn theorem and formal mathematical properties ... 
73. 21. La Quantum dots 
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I. INTRODUCTION 



Semi-relativistic Current-Density Functional Theory (CDFT) jsl, has 
become one on the standard tools for the calculation of electronic ground 
state (GS) properties in magnetic fields. Apart from the practical issue 
of finding an accurate and manageable energy functional, the basic fact 
of non-interacting V-representability (NIVR) of the exact electron density 
n(r) and the paramagnetic current density j p (r) is a prerequisite for the 
existence of a Kohn-Sham scheme. In a first step, we are investigating if 
both exact densities derived from the correlated two-particle state can be 
represented by the same set of one-particle orbitals 

Mr) = Rk(r) e <fe(r) (1) 
in the form (atomic units H = m = e = 1 are used throughout) 



occ occ 



n(r) = ^* fc (r) ^ fc (r) = 5>*(r) (2) 

k k 

1 occ occ 

j p ( r ) = 2; £ h (r) v ^ fc(r) " Mr) v ^ (r) ] = £^w ( 3 ) 

fc k 

where the modulus Rk(r) > and the phase Cfc( r ) °f the Kohn-Sham wave- 
functions (KS-WF) are real functions. If the KS-WFs are well defined, then 
the effective potentials V e -^ = (v 6 ^ and A e ^) could be obtained from the 
KS equations as exercised in Ref. |12j . In this paper we are only discussing 



if and when KS-WFs exist. Unlike in Density Functional Theory (DFT), 
NIVR in CDFT as defined above is neither guaranteed for nor restricted 
to GSs. 

For vanishing magnetic field, DFT applies, which rests on some unique 
mappings. If ^ describes a many-body GS-WF and n(r) and v ext (r) the 
corresponding density and external potential (modulo a constant), then of 

course ^ — > n(r) is unique, but also ^ — > v ext (r) , if infinitely high potential 
walls are excluded. Hohenberg and Kohn [5J proved that the mapping 

between ^ and n(r) is even one-to-one, and hence n(r) — > v ext (r) is also 
unique 

v ext (r) S- \& i n(r) (4) 

The situation of (Jlj) holds for both the interacting and the non-interacting 
cases. For a NIVR density n(r), the mapping CV' 1 for the non-interacting 
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case yields the (effective) KS potential of the interacting system. It has 
been shown (e.g. in [6(]) that not every mathematically well behaved den- 
sity n(r) is GS density to some potential v ext (r). Therefore, DFT has been 
based on functionals which are also defined for non-t>-representable den- 
sities. Nevertheless, a KS potential (derivative of the density functional) 
can only exist for NIVR densities. 

In the presence of a magnetic field and for (semi-relativistic) Current 
Density Functional Theory (CDFT), the generalization of V~ l for the 
ground state still exists, but, Vignale and Rasolt 0, 0] just presupposed 
the existence of the generalization of C [7| implying that NIVR and the 
existence of a KS scheme has not been proven. Capelle and Vignale |7|, on 
the other hand, have shown that there can be several external potentials 
yext w j 1 j c j 1 provide the same wave functions and densities 

Vf'(r) \ 

Vi xt (r) -»• AT(r) (5) 



where V ext (r) and J\f(r) represent both external potentials (v ext (r) and 
A ext (r)) and both densities (n(r) and j p (r)), respectively. Hence, C cannot 
exist anymore as a unique mapping. For Spin Density Functional Theory 
SDFT) the same problem was first pointed out by von Barth and Hedin 



| 111], and later analyzed in detail in [8|, |9j. For our model system the l.h.s. 
of (|I2) is obvious from the fact that the exact densities Af(r) are determined 
by the effective frequency 



u = yH 2 + {uj c /2f (6) 

alone and not by the external confinement frequency and the cyclotron 
frequency uj c = B/c independently (see Sect. II). In other words, all combi- 
nations of luq and cj c , which provide the same u, provide the same densities. 
This fact rules the existence of the mapping C out, but does neither prove 
nor rule out NIVR or the existence of a KS system. 



Wensauer and Rossler |12j used the scaling property of our quantum dot 
model system 

H(u , uj c ) = H(uj Q -> u, u c — > 0) + y L z , (7) 

where L z is the total orbital angular momentum, in order to apply the con- 
sequences from the Hohenberg-Kohn theorem to non-zero fields. Indeed, 
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([7j) means that the Hamiltonian for a non-zero magnetic field H(uo, cu c ) has 
the same eigen-functions as the Hamiltonian for zero magnetic field and 
the effective confinement frequency uj. Only the eigenvalues are shifted by 
(u c /2)Ml, where Ml is the total orbital angular momentum. The point is, 
however, that this does not mean that the GS densities J\f(r) for all u c (for 
fixed cjq) are NIVR. Instead, this conclusion applies only to those cases, 
where the corresponding reference state at zero uj c is a ground state. This 
are the states with Ml = 0, which are special insofar, as the paramagnetic 
current vanishes even for non-zero magnetic field. Their argument does 
not say anything about NIVR of the other cases. 

In the present paper we investigate, under which conditions the exact 
GS densities J\f(r) of our model system can be represented by KS orbitals. 
For making the conclusions gauge independent, we also matched the gauge 
independent vorticity 



instead of the paramagnetic current density j p (r). We will show that for a 
typical quantum dot model (two electrons in a parabolic confinement and a 
magnetic field) the exact GS densities are not generally NIVR. This is not a 
sophistry or subtleness, but the violation of NIVR is massive. This suffices 
to disprove the assumption of general NIVR and the general existence of 
a KS system for the GSs. On this background, all semi-relativistic CDFT 
calculations and functionals have to be considered with caution. On the 
other hand, this does not mean that all results are utterly wrong. However, 
this fact undermines the credibility of semi-relativistic CDFT as a tool for 
making forecasts. 

II. SPECIFICATION OF THE MODEL AND EXACT DENSITIES 
A. Model Hamiltonian 

We consider a two-dimensional two-electron system (with Coulomb in- 
teraction between the electrons) in a harmonic scalar potential v ext (r) = 
(1/2) Uq r 2 and a magnetic field perpendicular to the plane B = Be z rep- 
resented by the vector potential (in symmetric gauge) A ext (r) = (1/2) B x 
r = (1/2) B r e a . We introduced cylinder coordinates (r, a, z) with the 




(8) 
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cylinder axis perpendicular to the plane. The Hamiltonian reads 

^)V+>o^KuT— 77 (9) 




This is a widely used effective Hamiltonian model for a two-electron quan- 
tum dot. The interaction of the spins with the magnetic field H sp i n = 

3 

g* Si ■ B is omitted (by chosing g* = 0) for two reasons. First, in 

i=l 

semi-relativistic effective Hamiltonian theory g* is a material dependent 
parameter well below the vacuum value g = 2. Therefore, the limiting 
case of vanishing g* ought to be covered by an exact theory. Second, the 
Zeeman term would make the Hamiltonian spin dependent and would ne- 
cessitate a description of the system by the spin density n s (r) (instead 
of the total density n(r)) and the paramagnetic current density j p (r) Q. 
The introduction of the spin density in SDFT produces its own problems 
8, sl, 10] even for vanishing magnetic field, which we do not want to let 
interfere with the problems produced by the magnetic field. 



B. Exact solutions of the Schrodinger equation 



The Schrodinger equation with the Hamiltonian (J2j) can be solved 
not only by reduction to the numerical solution of an (ordinary) radial 
Schrodinger equation [E^] , but even analytically for a discrete, but infinite 
set of effective frequencies u {]], 0] . If we introduce relative and center of 
mass coordinate 



r = r 2 - i"! 



R=^(ri+r 2 ) 



the Hamiltonian (J9J) decouples exactly. 

1 



H = 2 H r + - H R + H spin 



(10) 



The Hamiltonian for the cm. motion agrees with the Hamiltonian of a 
non-interacting quasi particle 

2 



H 



R 



c 



1 

— ( 
2 



(12) 



and only the relative Hamiltonian contains the electron-electron interaction 

2 



H r = — 
2 



1 A 
p + -A 

c 



1 2 2 1 

+ 2^ r + Tr 



(13) 
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where we introduced rescaled parameters ujr = 2cjo, A# = 2A(R), u r = 
^cjo, A r = |A(r) (the index V and 'i?' refers to the relative and cm. 
coordinate systems, respectively). The decoupling of H allows the ansatz 

$ = «R)¥>(r)xK*2) (14) 

where x( s i> s 2) are the singlet or triplet spin eigen-functions. 
The eigen-functions of the cm. Hamiltonian ( I12j) have the form 

t = e -j= ^l = e R M (R) ; M = 0,±l,±2,... (15) 



where the polar coordinates of the cm. vector are denoted by (R, A) 
and the radial functions Um(R) and Rm{R) can be found in standard 
textbooks. 

With the following ansatz for the relative motion 



e ima Umir, 



^ = -= ; m = 0,±l,±2,... (16) 



the Schrodinger equation H r cp(r) = e r ip(r) gives rise to a radial 
Schrodinger equation for u(r) 

d 2 1 ( o 1\ 1 1 ~ 2 2 

-lj^ + 2^ r + 2rj 





>w(r) = e r u(r) (17) 



where the polar coordinates for the relative vector are denoted by (r, a), 
cu r = ^uj, ? r = e r — jiriLUc, and cu c = — . The solutions are subject to the 



normalization condition J dr\u(r)\ 2 = 1. The Pauli principle demands that 

o 

(because of the particle exchange symmetry of the spin eigen-functions) in 
the singlet and triplet state, the relative angular momentum m has to 
be even or odd, respectively. There is no constraint for the cm. angular 
momentum M following from Pauli principle. Because of the orthogonality 
of the coordinate transformation, the above described solutions are eigen- 
functions of the total orbital angular momentum with the eigenvalue Mi = 
M + m. 

FigH] shows that the modulus of the orbital angular momentum of the 
ground state grows stepwise with increasing magnetic field. This implies 
that the spin state oscillates between singlet and triplet [14]. (Quenching 
of the singlet state for higher magnetic fields due to a Zeeman term is not 
included in our model.) States with cm. excitations are not included in 
the figure, because they are never ground states. 
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FIG. 1: (color online) Total energy for fixed confinement frequency Uq = 1 versus cyclotron 
frequency lo c (i.e. magnetic field). The cm. system is always in the ground state with M = 0. 
The relative angular momentum m is varied. The vertical lines show where the total orbital 
angular momentum Ml = M + m of the ground state changes. S and T indicates whether the 
ground state is singlet or triplet. Thick lines indicate states which can be NIVR. 



C. Exact densities 

With ( |T5| ) and (JTBT) . we obtain for the total density 

./\|2 



rc(r) = 2 J dr' |$(r,r ; )| : 



(18) 



the general expression 
1 



n(r) = 



2tt 2 



' da / dr 1 

o Jo 



1 

4' 



Rm ( \ / r 2 + -r' 2 + rr'cosa 



u m {r') 



(19) 

Because we are interested in the ground state only, we can safely use the 
cm. state for M = 0: Ro(R) = 2Vcu exp(—cuR 2 ) which allows to do one 
integration analytically leaving us with 



n{r) = — e -^ r2 H dr' e^ 2 ^ I (2urr') 

7T Jo 



u m (r') 



(20) 



where I n (x) are the modified Bessel functions. 
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The general expression for the paramagnetic current density 



f(r) = -i dr' $*(r,r')V$(r,r') - $(r, r')V$*(r, r'} 



(21) 



is somewhat complicated. Therefore we give here only the formula for 
M = 



f(r; 



4cj 



e a m 



-2wr 2 



dr' e-^ 2 > r 



h(2urr') 



7T 



Um\r 



e« fir) 



(22) 

As to be expected, the paramagnetic current density points in azimuthal 
direction e a , and the scalar j p (r) depends only on the distance r from the 
center and not from the azimuthal angle . 

Although both formulas ( 1201 ) and ( 1221 ) rely on the functions u m (r), which 
are solutions of (JT71). the analytical behavior for r — » can be expressed in 



S 



terms of two positive definite integrals 



A, = 
A 2 = 



dre-M** 1 * u m (r) 



dr r 2 e-^' 2 > 2 



u m {r) 



(23) 
(24) 



After power series expansion of I n (x), we obtain 

4uj 



n (r 



f(r) 



TT 

m 



A + A 2 uj 2 r 2 + 



4u 2 



TT 



1 

2 



A + ^A 2 uj 2 r 2 + 



(25) 
(26) 



For the origin this means that n(0) = 4 to Aq/tt is always finite and j p (0) = 
always vanishes. On the other hand, the derivative of the density at 
the origin 4^(0) = vanishes, but of the paramagnetic current density 



-Jr(0) = m(4 to 2 1 V) Aq is finite, unless m = 0. Besides, there is a relation 
which does not rely on the radial functions at all 



dr 



■(0) = m uj n(0) 



(27) 



The exact vorticity, which has the form 'y(r) = e z 7(7"), reads in this limit 

'■ 2 ' (28) 



m 2 cu [ 1 — uj" r" + 
^4o 



As will be seen below, the limit r — ► is decisive for our proof. The overall 
behavior of the densities is shown in FigJ21 

III. DETERMINE KOHN-SHAM ORBITALS FROM EXACT DENSITY AND 
PARAMAGNETIC CURRENT DENSITY 



Our consideration is in a sense complementary to the approach in |12 



where the triplet GS with Mi = — 1 was considered. We are going to point 
out that this is the only triplet GS for which can be NIVR. 

A. Singlet state: 

In the KS system we assume that both electrons occupy the same or- 
bital state of the form (p(r) = R(r) with different spins. We do not 
presuppose that this is an eigenstate of the orbital angular momentum, 
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because there is no theorem which proves that orbital angular momentum 
is conserved in the KS system, if it is conserved in the real system. Now we 
demand that both densities (n, } p ) in the KS and in the real system agree. 



n KS (r) = 2 [fl(r)] 2 = n exact (r) (29) 



f KS (v) = 2 [R(v)f VC(r) A f exact (r) e a (30) 
Equation ( 1291 ) shows that the real part depends only on r and reads 



R(r) = J ^ n exact (r) (31) 
Inserting the gradient in polar coordinates 

V((r,a) = e r — + e a -— (32) 
into ( 1301 ) and using ( jZSD provides two equations for £(r, a): 

%^A _ f = ((a) (33) 

= T = /lexactir) Vexact{r) = COUSt. (34) 

eta n exac t{r) 

If a function of a has to agree with a function of r (left equation of (1341)). 
then both sides have to be a constant. Consequently, /or NIVR the function 
Hexactir) has to be constant. 

Before investigating this issue further, we introduce the vorticity (jSj) of 
the exact system. By using the special form of the curl of a vector field v, 
which points in e a direction and the modulus of which depends only on r, 

V x v(r) e a = -— ( r v(r) )e z (35) 
r dr \ J 

we see that the vorticity of the exact solutions has the general form 

lexacti?) = ~ ( r^^) e z = Jexact (r) e z (36) 
rdr \ n exact (r) J 

where e z is the unit vector perpendicular to the 2D system, and from ( 1341 ) 
it follows 

lexact(r) = -— Hexact(r) = (37) 

r dr 

This means that for NIVR the exact vorticity 7 exac t(r) has to vanish. 



10 



The pure fact of violation of NIRV (without showing the quantitative 
extent in r-space) can already be seen analytically in the limit for small 
r given in ( 1251) . which provides 7 e xaci(0) = 2 to Ml. This value can be 
considered as a simple qualitative measure for the degree of violation. The 
numerical curves in Figs. El and [U show that the violation of NIVR in 
an extented small-r-region is massive. There is only one singlet state, 
where NIVR is trivially guaranteed for all r. This is Ml = 0, where 
the paramagnetic current and consequently fiexacti^) and jexacti^) vanish 
exactly. (As a side product we observe that fi eX act{r) converges for large 
r to a constant, which agrees with the (total) orbital angular momentum 
Ml of the exact state.) 

Next we are investigating the vorticity of the KS system. From ((HI) 
and ((2211313]) it follows that the vorticity of one general KS state vanishes 
automatically 



irrespective of its special form. If two electrons occupy the same orbital 
state then the density and paramagnetic current density have to be multi- 
plied by the factor of 2 and the vorticity of the whole KS system vanishes 
automatically. 

One could object that the above described procedure is not gauge invari- 
ant, because j p (r) is gauge dependent. If we replace ( 15U1 ) by equating the 
gauge invariant vorticities of both systems and consider that the vorticity 
of the (singlet state in the) KS system vanishes 



it is obvious that the exact vorticity has to vanish, what agrees with the 
result of the former approach. However, in order to be absolutely correct 
we have to consider the following subtlety. Observe that the canonical 
orbital angular momentum ^ k x is gauge dependent and it is only 
conserved in the symmetric gauge (see Sect.IIA), which we used in out 
exact solutions. Therefore, we should characterize the only singlet state, 
which can be NIVR, in the following way: it is that state, which has Ml = 
in the symmetric gauge. This statement covers all gauges. 




(38) 




(39) 
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FIG. 3: l^exacti^) (in a.u.) for loq — 1 and typical cyclotron frequencies u> c , where the state with 
negative Ml is the ground state (see Fig(T]). The sign of ^ e xact{f) agrees with the sign of Ml- 
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FIG. 4: Vorticities (in a.u.) for the same external fields as Figj3j The sign of 7 eX act(^) agrees 
with the sign of Ml- 



B. Triplet state: 



In a first step of sophistication we assume that the KS system is circular 
and the KS states are eigen-functions of the orbital angular momentum 
<£>fc(r) = e imka Rk(r) . (Do not mix up the relative angular momentum m in 
in Sect. 2 and the one-particle angular momentum in this section, which 
are both denoted by 'm'.) Therefore, we investigate, whether the exact 
system (which is circular and conserves total orbital angular momentum) 
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can be replaced by a circular KS system with the same total angular mo- 
mentum. This assumption is common practice in numerical self-consistent 



cLS cl 



CDFT calculations and allows us to obtain the results in Ref. [12 
special case. Equating both densities (n, } p ) provides the following set of 
equations 



n KS (r) = [R^r)} 2 + [R 2 (r)} 2 = n exact (r) (40) 

Iksv ) - m i r- m 2 J exact {r) (41) 

Now we investigate the limit for small r. In the Appendix it is shown that 
for the special case of a circular KS system Rk{r) — > Ck r' TOfc L Equations 
(E5j) and (ESD can be written as n exact (r) —> no + n 2 r 2 and f exact (r) — > 
jo Ml r. Then (SDj) and ( l41l ) read in the limit r — > for the essential 
lowest terms 

c 2 r 2K| + c 2 r 2|m 2 | =n() + n2r 2 ^ 

c\ mi r 2 ^- 1 + c 2 2 m 2 r 2 '" 12 '" 1 = j M L r (43) 

( I42D can be fulfilled only if one of the angular momenta vanishes and the 
modulus of the other one is unity, say m\ = and \m 2 \ = 1. This choice 
satisfies also (1431). Consequently, only the exact states with \Ml\ = \m\ + 
7?t-2 I = 1 can &e NIVR by a circular KS system, whereby only the state 
with Ml = — 1 can be the ground state (see FigJT]). In this case we obtain 
from ( ]4"0j ) and ( 14~TT ) explicit formulae for the radial parts 

[^(r)] 2 = -r^ act (r) (44) 

[i?i(r)] 2 = n exaci (r) - [R 2 (r)} 2 (45) 



agreeing with Ref. fl2| 



In a second step, we abandon the constraint for circular symmetry of 
the KS system, because it cannot be taken for granted. Additionally, we 
equate the gauge invariant vorticities instead of the gauge-dependent para- 
magnetic current densities of the exact and the KS system. The analytical 
properties of the KS functions for r — ► in this general case are discussed 
in Appendix A. The gauge constants rag can be neglected here because 
they cancel in gauge invariant quantities like n(r) and 'y(r) anyway. If we 
substitute Rk(r) — > Ck r mk in both KS densities ( l4"0j ) and (14T|). and insert 
both densities into the the vorticity of the KS system 




j ks (t) = V x ea = [ r ^± Gz = 7 ( r ) e , (46) 
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we obtain 

2c\ c\ (mi - m 2 ) 2 r 2 ^ +m ^ 

lKS\f) ~^ ^2 r 2mi c 2 r 2m 2 j2 V 

The lowest-power term of the exact vorticity ( 1251 ) can be written as 

7™t(r) 70 M L (48) 
Consequently, the two constraints from equating n(r) and j(r) read 

c 2 r 2m! + c 2 r 2m 2 = ^ + ^ r 2 ^ 

2 c\ c\ (mi - m 2 ) 2 r 2 *™^" 1 ) . _ 

? ■ 2T2 = ^° Mi 50 ) 

whereby in the denominator of the second constraint the first constraint 
has been used. Although ( 150 ) and ( 1431 ) differ, the conclusion from ( 14211431 ) 
and ( 1491150 ) are very similar. ( 140150 ) can be fulfilled only if the small-r- 
exponents of the KS states are and 1 . The same conclusion can be drawn 
from equating the density and the paramagnetic current density for KS 
systems of arbitrary symmetry. 
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IV. SUMMARY AND CONCLUSIONS 



We used the exact solutions of the Schrodinger equation of a two- 
dimensional model system for the investigation of NIVR. Qualitative con- 
clusions can be drawn already on the basis of analytical results for the 
power expansion for r — » of the densities J\f{r). 

For singlet states {Ml = even) there is a simple and straightforward 
proof without any assumptions that only the state with vanishing total 
orbital angular momentum Ml = can be NIVR. This state is the ground 
state for small magnetic fields (see FigJTJ) and it does not produce a para- 
magnetic current. 

For triplet states {Ml = odd) a fully satisfactory statement can be made 
only under the assumption, that the KS system has the same circular 
symmetry as the real system and consequently conserves orbital angular 
momentum. Then only the exact states with \Ml\ = 1 are NIVR by KS 
states with angular momenta ni\ = and \ni2\ = 1- 

If we allow KS systems with arbitrary symmetry, we can show that NIVR 
allows only KS states with small-r-exponents (see Appendix) of and 1. 
However, we cannot say which exact states can or cannot be represented 
by symmetry-broken KS systems. In other words, we know that the exact 
states with \Ml\ = 1 can be represented by circular KS systems, and 
that all other states cannot be represented by circular KS systems. This 
statement is valuable despite its limitations, because the assumption of 
circular symmetry in self-consistent calculations is common practice. 

At the end we want to discuss the connection between NIVR and the 
property of being a ground state. The singlet states belonging to the full 
black line {Ml = 0) in FigQ] can be NIVR, no matter whether they are 
ground states for the given coq and uj c or not. All other states with even 
Ml cannot, even if they are ground states. The full blue line {Ml = — 1), 
which is the ground state in the second o; c -region in Fig HI can be NIVR 
everywhere. Even the states belonging to the broken blue line {Ml = 4-1) 
are NIVR, although they are never the ground state. All other states with 
odd Ml are not NIVR by a circular KS system. Consequently, all those 
states at non-zero B can be NIVR, which are continuously connected to 
states at B = 0, which are the ground states for a given spin state at 
B = 0, no matter if they are the ground states at non-zero B or not. The 



last statement sounds similar to the scaling property in Ref. fl2j , but it is 
not identical (see Introduction). 
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APPENDIX A: GENERAL KOHN-SHAM WAVE FUNCTION FOR r 
AND GAUGE DEPENDENCE 



We are going to show that (apart from a gauge-invariant normalization 
factor) any KS wave function in the limit r — > has the form 



^Pks\T,ol j 



(Al) 



where the small-r- exponent m > is a gauge-invariant integer, but the 
integer mg depends on the gauge. If the canonical angular momentum 
l z = (l/i)d/da is conserved, then fh + mg is the canonical angular momen- 
tum m of the KS state. Otherwise it is just a state and gauge-dependent 



integer. The limiting behavior of the density n — > r 
as expected, but the paramagnetic current density j p 



2m 



is gauge- invariant, 



,777, 



mg) r 



2m- 1 



depends on the gauge. Even if the canonical angular momentum is con- 
served, the exponent determining the radial WF (and the density) for small 
r and the canonical angular momentum should not be mixed up. 
The KS Hamiltonian has the general form 



Hrs = 



1 



-V + ^A(r) 

i c 



+ V(r) 



(A2) 



The concrete form of the effective potentials in this equation and their 
connection to the XC-energy functional can be found in |4j • Now we apply a 
gauge transformation with x — —crag a which changes the vector potential 
as follows 

1 

A -> A' = A + (Vx) = A - cm g -e a (A3) 

r 



This gauge transformation is equivalent to introducing a flux line with mg 
flux quanta at the center. After rearranging terms and writing the Laplace 
operator in polar coordinates we obtain 



H' = 



i r -i/2^ 



T l/2 + 



1 



1 A 

H mg — e 

c r 



lei 



ld_ 

i da 



mg 



1 



+ — A • V + —jA z + V (A4) 
ci Zc l 



Because in our model all external potentials are continuous, it is natural 
to assume that so is V and A. (Actually we have only to assume that V 
diverges weaker than 1/r 2 and A weaker than 1/r.) Then the eigenfunction 
in the limit r — > is defined by the first line in (IA4j ) . which is independent 
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of the specific form of V and A. With the ansatz 

u(r) 



ip = 



r 



1/2 



w(a) 



and after multiplication with r 2 [l/u{r)w{a)]r l l 2 the variables in the 
Schrodinger equation H' ip' = e if' can be separated providing 



u(r) 



1 2 d 2 
—r 

2 Or 2 



r 2 e 



u(r) 



1 1 



2 w(a) 




mg + 



w (a, 



(A5) 



This equation can only be fulfilled if both sides are constant. Using the 
technique of separation of variables, we can easily show that (ED) is the 
solution, whereby the term with r 2 e on the l.h.s. can be neglected because 
we are interested in the limit r — > 0. 
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